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The constituent parts of a quantum computer are inherently vulnerable to errors. To this end 
we have developed quantum error-correcting codes to protect quantum information from noise. 
However, discovering codes that are capable of a universal set of computational operations with the 
minimal cost in quantum resources remains an important and ongoing challenge. One proposal of 
significant recent interest is the gauge color code. Notably, this code may offer a reduced resource cost 
over other well-studied fault-tolerant architectures using a new method, known as gauge fixing, for 
performing the non-Clifford logical operations that are essential for universal quantum computation. 
Here we examine the gauge color code when it is subject to noise. Specifically we make use of 
single-shot error correction to develop a simple decoding algorithm for the gauge color code, and we 
numerically analyse its performance. Remarkably, we find threshold error rates comparable to those 
of other leading proposals. Our results thus provide encouraging preliminary data of a comparative 
study between the gauge color code and other promising computational architectures. 


INTRODUCTION 

Scalable quantum technologies require the ability to 
maintain and manipulate coherent quantum states over 
an arbitrarily long period of time. It is problematic 
then that the small quantum systems that we might use 
to realize such technologies decohere rapidly due to un¬ 
avoidable interactions with the environment. To resolve 
this issue we have discovered quantum error-correcting 
codes mu which make use of a redundancy of physical 
qubits to maintain encoded quantum states with arbi¬ 
trarily high fidelity over an indefinite period. 

Ideally, we will design a fault-tolerant quantum com¬ 
puter that requires as few physical qubits as possible to 
minimize the resource cost of a quantum processor, and 
indeed, the cost in resources of a computational archi¬ 
tecture is very sensitive to the choice of quantum error- 
correcting code used by a fault-tolerant scheme. It is 
therefore of great interest to analyse different quantum 
error-correction proposals to compare and contrast their 
resource demands. 

Color codes GMZ] are a family of topological quan¬ 
tum error-correcting codes [SHm with impressive versa¬ 
tility pjhh] for performing fault-tolerant logic gates da- 
ini. This is an important consideration as we search for 
schemes that realise fault-tolerant quantum computation 
with a low cost in quantum resources. In particular, a 
fault-tolerant quantum computer must be able to per¬ 
form a non-Clifford operation, such as the 7r/8-gate, to 
realise universal quantum computation. 

In general, performing Non-Clifford gates can present 
a considerable resource cost over the duration of a quan¬ 
tum computation. As such, the resource cost of real¬ 
ising scalable quantum computation is sensitive to the 
method a fault-tolerant computational scheme uses to 


realise non-Clifford gates. To this end the gauge color 
code [3 [H3-ED3 has attracted significant recent interest 
because, notably, this three-dimensional quantum error- 
correcting code can achieve a universal gate set via gauge 
fixing pH! [221 . 

In contrast to the gauge-color code, surface code 
quantum computation, a leading approach towards low- 
resource quantum computation [3 13 US], makes use of 
magic state-distillation [24] to perform 7r/8-gates. Magic 
state distillation can be achieved with 0(L 3 ) space-time 
resource cost [23J I25H25] . Similarly, the gauge color code 
performs 7r/8-gates via gauge fixing in constant time [T91 . 
and as such has an equivalent scaling in space-time re¬ 
source cost as the surface code since the gauge color 
code requires 0{L 3 ) physical qubits. However, given 
that gauge fixing requires no additional offline quantum 
resources to perform a non-Clifford rotation, the gauge 
color code may reduce the quantum resources that are 
necessary for fault-tolerant quantum computation by a 
constant fraction. 

It is also noteworthy that the gauge color code is local 
only in three dimensions and as such, unlike the sur¬ 
face code, cannot be realised using a two-dimensional 
array of locally interacting qubits. Instead, the gauge 
color code may be an attractive model for non-local 
quantum-computational architectures such as networked 
schemes [291133] , 

Given the significant qualitative differences between 
the gauge color code and the surface code, it is interesting 
to perform a comparative analysis of these two propos¬ 
als. In this Manuscript we investigate error correction 
with the gauge color code. 

Dealing with errors that continually occur on physi¬ 
cal qubits is particularly difficult in the realistic setting 
where syndrome measurements can fail and return false 
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readings [5]. Attempting to correct errors using inaccu¬ 
rate syndrome information will introduce new physical 
errors to the code. However, given enough syndrome in¬ 
formation, we can distinguish measurement errors from 
physical errors with enough confidence that the errors we 
introduce are few, and can be identified at a later round 
of error correction |T5] . In the case of the toric code [9], 
we accumulate sufficient error data by performing mul¬ 
tiple rounds of syndrome measurements. Surprisingly, 
the structure of the gauge color code enables the acquisi¬ 
tion of fault-tolerant syndrome data using only one round 
of local measurements nu. This capability is known as 
single-shot error correction. 

Here, we obtain a noise threshold for the gauge color 
code using a phenomenological noise model where both 
physical errors and measurement faults occur at rate p. 
We develop a single-shot decoder to identify the sustain¬ 
able operating conditions of the code, i.e. the noise rate 
below which information can be maintained arbitrarily 
well, even after many cycles of error correction. We es¬ 
timate a sustainable error rate of p sus ~ 0.31% using an 
efficient clustering decoding scheme (34H391 that runs in 
time 0(L 6 log L), [35] . where the distance of the code is 
d = L + 2, [7|. Remarkably, the threshold we obtain falls 
within an order or magnitude of the optimal threshold for 
the toric code under the same error model, ~ 2.9%, m- 
Furthermore, we also use our decoder to estimate how the 
logical failure rate of the gauge color code scales below 
the threshold error rate by fitting to a heuristic scaling 
hypothesis. 

RESULTS 

The gauge color code 

The gauge color code is a subsystem code m speci¬ 
fied by its gauge group, Q. From the center of the gauge 
group, Z(Q), we obtain the stabilizer group for the code, 
S = Z(G) n G, and its logical operators, C = Z{G)\G- 
Elements of the stabilizer group, S G S, satisfy the prop¬ 
erty that S\ip) = | ip) for all codewords of the code 

The code is defined on a three-dimensional four-valent 
lattice of linear size L with qubits on its vertices mm- 
The lattice must also be four-colorable, i.e. each cell of 
the lattice can be given a color, red, green, yellow, or 
blue, denoted by elements of the set C = {r,g,y, b}, 
such that no two adjacent cells are of the same color. 
The lattice we consider is shown in Fig. [IJa), where the 
cells are the solid colored objects in the Figure. 

The cells of the lattice define stabilizer generators for 
the code, while the faces of cells define its gauge gen¬ 
erators. The gauge operators, otherwise known as face 
operators, are measured to infer the values of stabilizer 
operators. The face operators have weight four, six and 
eight, where the weight-eight face operators lie on the 



FIG. 1. (a) In the primal picture qubits sit on the vertices 
of a four-valent lattice. The three-dimensional cells of the 
lattice are four colorable, i.e., every cell can be assigned one 
of four colors such that it touches no other cell of the same 
color, (b) The gauge color code of linear size L = 5 drawn 
in the dual picture. Qubits lie on simplices of the lattice. 
The faces of the tetrahedra in the figure are given one of 
four colors such that no two faces of a given tetrahedron have 
the same color. The tetrahedra are stacked such that faces 
that touch always have the same color. For the gauge color 
code to encode a single logical qubit, the lattice must have 
four distinct, uniformly colored boundaries, as is shown in 
the Figure. We give more details on the lattice construction 
in the dual picture in Supplementary Note 1 together with 
Supplementary Figures 1-6. 

boundary of the lattice. More specifically, for each cell 
c there are two stabilizer generators, S* = IljeV(c) 
and Sf? = Jljev(c) where V(c) are the qubits on the 
boundary vertices of cell c, and X 3 and Z 3 are Pauli- 
X and Pauli-Z operators acting on vertex j. For each 
face / there are two face operators, Gj = IIjev(/) 
and Gj = IIjev(/) where V(f) are the vertices on the 
boundary of face /. We will call the outcome of a face 
operator measurement a face outcome. 

Given suitable boundary conditions emu, the code 
encodes one qubit, whose logical operators are X = 
J| j6 g Xj and Z = IljeC ■ where Q is the set of physi¬ 
cal qubits of the code. A lattice with correct boundaries 
is conveniently represented on a dual lattice of tetrahe¬ 
dra using the convention given in Ref. [Tj. The lattice 
we consider is illustrated in Fig. |T](b) , where we discuss 
its construction in detail in the Supplementary Notes 1 
and 2, together with Supplementary Figures 1-7. Impor¬ 
tantly, we require that the lattice has four boundaries, 
distinguished by colors from set C. 

Single-shot error correction 

A quantum error-correcting code is designed to iden¬ 
tify and correct errors. Due to the symmetry of the gauge 
group, it suffices here to consider only bit-flip, i.e. Pauli- 
X errors. We consider a phenomenological error model 
consisting of physical errors and measurement errors. A 
physical error is a Pauli-X error on a qubit, while a mea¬ 
surement error returns the opposite outcome of the cor- 
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rect reading. Errors will be identically and independently 
generated with the same probability p. 

In a stabilizer code, errors are identified by stabilizer 
measurements that return eigenvalue —1, which we call 
stabilizer defects. We use the stabilizer syndrome, a list 
of the positions of stabilizer defects, to predict the inci¬ 
dent error. In the gauge color code we do not measure 
stabilizer operators directly, but instead infer their val¬ 
ues by measuring face operators, which is possible due to 
the fact that S C Q. 

In addition to using face outcomes to infer stabilizer 
eigenvalues, we can also exploit the local constraints in 
Q of the gauge color code to detect and account for mea¬ 
surement errors. Remarkably, measurement errors can 
be detected reliably by measuring each face operator only 
once, so called single-shot error correction il9j. 

The local constraints stem from the structure of the 
code. We first observe that the faces of a cell are neces¬ 
sarily three colorable, as shown in Fig. |2ja). It follows 
from the three colorability of the faces that the product of 
the gauge operators Gj of any of the differently colored 
subsets of faces of cell c recover the stabilizer operator 
. By measuring all the faces of the lattice, we re¬ 
dundantly recover each stabilizer eigenvalue three times 
where the three outcomes of a given cell are constrained 
to agree. In Fig. [2jb) we show an example of a gauge 
measurement configuration where the outcomes are reli¬ 
able. Following this observation, we can use violations of 
the local constraints about a cell to indicate the positions 
of measurement errors. In Fig. [2jc) we show a syndrome 
where the product ofthe face outcomes of different sub¬ 
sets of of two cells, colored with thick green edges, do not 
agree. 

Fault-tolerant decoding with the gauge color code pro¬ 
ceeds in two stages. The first stage, syndrome estima- 



FIG. 2. (a) The faces of the cells of the color code are three 
colorable, (b) An example of a set of face measurement out¬ 
comes where measurements are reliable. Face operators that 
return value -1 are colored, otherwise they are left transpar¬ 
ent. The stabilizer at the cell with thick red edges contains a 
stabilizer defect. This cell has one face operator of each color 
returning a -1 outcome. All other cells have an even parity 
of -1 face outcomes over their colored subsets, indicating that 
these cells contain no stabilizer defect, (c) A gauge syndrome 
where colored subsets of face measurements about the green 
cells do not agree, thus indicating measurement errors. 


tion, uses face outcomes that may be unreliable to esti¬ 
mate the locations of stabilizer defects. The second stage, 
stabilizer decoding, takes the estimated stabilizer defect 
locations and predicts a correction operator to reverse 
physical errors. 

The topological nature of the gauge color code is 
such that its stabilizer defects satisfy conservation laws 
that enable us to employ well-studied decoding algo¬ 
rithms mmm to complete stabilizer decoding. The 
syndrome identified by stabilizer measurements is stud¬ 
ied extensively in Refs. [4j |42]. Stabilizers that return 
— 1 outcomes occur in pairs at the endpoints of strings 
of errors, whose endpoints take the color of its terminal 
cells, as shown in Fig. [3^a). As such, one can also regard 
an error string as carrying the color of the stabilizer de¬ 
fects at its endpoints. Error strings on the gauge color 
code can also branch into three strings of conjugate col¬ 
ors, thus creating one syndrome of each color, as shown 
in Fig. §b). Finally, error strings of a given color can 
terminate at the boundary of their respective color, as 
shown in Fig. [3^c) . In our simulation we decode sta¬ 
bilizer defects by adapting a clustering decoder [34U39] 
where clusters grow linearly. 

We concentrate now on syndrome estimation. Syn¬ 
drome estimation uses a gauge syndrome, a list of gauge 
defects, to estimate a stabilizer sydrome. A lattice cell 
can contain as many as three gauge defects. Gauge de¬ 
fects are distinguished by a color pair uv, with u, v S C, 
such that u ^ v and uv = vu. The color pair of a gauge 
defect relates to the coloring of the lattice faces. A face is 
given the color pair opposite to the colors of its adjacent 
cells, i.e. the face shared by two cells with colors r and g 
is colored yb. A cell c contains a uv gauge defect if the 
product of all the uv face outcomes bounding c is — 1. 
Following this definition, a stabilizer defect is equivalent 
to three distinct gauge defects in a common cell. 

Studying the gauge syndrome enables the identification 
of measurement errors. We consider face /, colored uv, 
that is adjacent to cell c. In the noiseless measurement 
case, where c contains no stabilizer defect, by definition, 
cell c should contain no gauge defects. However, if face / 
returns an incorrect outcome, we identify a uv gauge de¬ 
fect at c. Conversely, if cell c contains a stabilizer defect 
in the ideal case, and / returns an incorrect outcome, no 
uv gauge defect will appear in c. With these examples 
we see that cells that contain either one or two gauge 
defects indicate incorrect face outcomes. 

An incorrect face outcome affects gauge defects in both 
of its adjacent cells. In general, incorrect face outcomes of 
color uv form error strings on the dual lattice, whose end 
points are uv gauge defects, where individual incorrect 
face outcomes are segments of the string. Error strings 
of incorrect uv face outcomes changes the parity of uv 
gauge defects at both of its terminal cells. 

We require an algorithm that can use gauge syndrome 
data to predict a likely measurement error configuration, 
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and thus estimate the stabilizer syndrome. We adapt the 
clustering decoder [841488] for this purpose. The decoder 
combines nearby defects into clusters that can be con¬ 
tained within a small box. Clusters increase linearly in 
size to contain other nearby defects until they contain a 
set of gauge defects that can be caused by a measure¬ 
ment error contained within the box. Once clustering 
is completed, a correction supported inside the boxes is 
returned. 

We briefly elaborate on correctable configurations of 
gauge defects. Pairs of uv gauge defects are caused by 
strings of incorrect uv face outcomes, and therefore form 
correctable configurations, as shown in Fig. [3^d). As an 
example, the gauge syndrome in Fig. id) depicts the 
gauge defects and measurement error string shown in 
Fig. [2jc) where the measurement errors have occurred 
on the faces that returned —1 measurement outcomes. 
Triplets of gauge defects, colored uv, vw and uw, also 
form correctable configurations. Error strings that cause 
triplets of correctable gauge defects branch at a cell, and 
thus indicate a stabilizer defect, shown in Fig. He). The 
stabilizer defect where the error string branches lies at 
a cell colored x ^ u,v,w. Gauge defects can also arise 
due to incorrect face outcomes on the lattice boundary. 
Specifically, the boundary colored w contains faces of 
color uv where u,v/w. With this coloring we can find 
correctable configurations of single uv gauge defects, to¬ 
gether with a boundary of color w, see Fig.[3]T). In gen¬ 
eral, a cluster can contain many correctable pairs and 
triplets of gauge defects. 

As we have mentioned, correctable clusters of gauge 
defects can give rise to stabilizer defects. It is important 
to note that the error-correction procedure is sensitive 
to the positions of stabilizer defects within a correctable 
cluster, as discrepancies in their positions later affect the 
performance of the stabilizer decoding algorithm. As 
such we must place stabilizer defects carefully. For cases 
where a correctable cluster of gauge defects returns sta¬ 
bilizer defects, we assign their positions such that they 
lie at the mean position of all the gauge defects within 
the correctable cluster, at the nearest cell of the appro¬ 
priate color. Once syndrome estimation is complete, the 
predicted stabilizer syndrome is passed to the stabilizer 
decoder, and a correction operator is evaluated. 

We remark that gauge defects can be incorrectly ana¬ 
lyzed during syndrome estimation. In which case, mea¬ 
surement errors sometimes masquerade as stabilizer de¬ 
fects, and sometimes stabilizer defects can be misplaced. 
We will then attempt to decode the incorrect stabilizer 
syndrome and mistakenly introduce errors to the code. 
In general, any error-correction scheme that takes noisy 
measurement data will introduce residual physical errors 
to a code. These errors can be corrected in the future, 
provided the remaining noise is of a form that a decoder 
can correct. In general however, one must worry that 
large correlated errors can be introduced that adversar- 



FIG. 3. (a) Two blue syndromes indicate the end points of a 
string of errors connecting the two points, (b) A string error 
branches to create one syndrome of each color, (c) A green 
string error can terminate at a green boundary, thus gener¬ 
ating only a single green syndrome, (d) A pair of yb gauge 
defects, shown by the vertices, can be caused by a string of 
incorrect yb face outcomes on the dual lattice. The displayed 
gauge syndrome is equivalent to that shown in Fig.[2[c) where 
the measurement errors have occurred on the faces that re¬ 
turned — 1 outcomes, (e) Three gauge defects, colored gy, gb 
and yb, can be caused by an error string that branches at a 
red cell. The branching point indicates a stabilizer defect at 
a red cell, (f) Strings of incorrect face outcomes of colors rg 
and gb terminate at a yellow boundary. 

ially corrupt encoded information [451450] . Such errors 
may occur in the gauge color code if, for instance, we 
mistakenly predict two stabilizer defects of the same color 
that are separated by a large distance. We give an exam¬ 
ple of a mechanism that might cause a correlated error 
during syndrome estimation with the gauge color code in 
the Supplementary Note 3, together with Supplementary 
Figure 8. 

A special property of the gauge color code is that mea¬ 
surement errors, followed by syndrome estimation, will 
only introduce false defects in locally correctable config¬ 
urations. Therefore, residual errors remain local to the 
measurement error. Moreover, the code is such that the 
probability of obtaining configurations of face outcomes 
that correspond to faux stabilizer defects decays expo¬ 
nentially with the separation of their cells. This is be¬ 
cause the number of measurement errors that must occur 
to produce a pair of false stabilizer defects is extensive 
with their cell separation. To this end, the errors intro¬ 
duced from incorrect measurements are local to the mea¬ 
surement error and typically small. This property, coined 
‘confinement’ in Ref. m, is essential for fault-tolerant 
error correction. Most known codes achieve confinement 
by performing syndrome measurements many times. We 
give numerical evidence showing that our error-correction 
protocol confines errors in the following Subsection. 

The Simulation 

We simulate fault-tolerant error correction with en¬ 
coded states | ipj) of linear size L where j indicates the 
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FIG. 4. We show logical error rates, Pfaii, as a function of 
physical error rate, p, for system sizes L = 17, 23, 29 and 
35 shown in blue, yellow, green and red, respectively, as is 
marked in the legend, where we collect data after N = 8 
rounds of error correction during which measurements are 
performed unreliably. The error bars show the standard error 
of the mean given by the expression A p = \/{l — P{a,i\)Pui\/p 
where p is the number of Monte Carlo samples we collect. 
The data used to determine the threshold error rate is shown 
in the inset, where we determine the threshold using the fit¬ 
ting described in the Methods Section. The fitting is also 
plotted in the inset. In the main Figure the solid lines show 
the fitted expression, Eqn. 0. to demonstrate the agreement 
of our scaling hypothesis with numerically evaluated logical 
error rates where p < pth- We remark that the fitting is made 
using data for all values of TV, and not only the data shown 
in this plot, as we explain in the Methods Section. 

number of error-correction cycles that have been per¬ 
formed, and where |^o) is a codeword. We seek to find a 
correction operator C such that CE £ Q where E is the 
noise incident to j'i/’o) after TV error-correction cycles. 

To maintain the encoded information over long du¬ 
rations, we repeatedly apply error-correction cycles to 
keep the physical noise sufficiently benign. After a short 
period, the state \ipj-i) will accumulate physical noise 
Ej ( p ) with error rate p. To correct the noise, we first es¬ 
timate a stabilizer syndrome, s,, using gauge syndrome 
data with the syndrome estimation algorithm M q . where 
measurement outcomes are incorrect with probability 
q = p. Specifically, we have that Sj = M p (Ej(p)\ipj_i)). 
We then use the stabilizer decoding algorithm D to pre¬ 
dict a suitable correction operator Cj = D(sj ), such that 
we obtain 

= (!) 

It is important to note is that | ipj) is not necessarily in 
the code subspace. For q > 0, stabilizer syndromes will in 
general be incorrectly estimated, and thus the correction 
operator Cj introduces some new errors to the code. 

We require that after TV error-correction cycles we can 
estimate error E of state \ipt?) = E\i/jq) to perform a 
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FIG. 5. Threshold error rates, pth, are calculated with sys¬ 
tem sizes L = 23, 29, 35 after TV error-correction cycles using 
p ~ 10 4 Monte Carlo samples. Error bars show the standard 
error of the mean which are determined using the NonLin- 
earModelFit function in Mathematica. The solid blue line 
shows the fitting given in Eqn. 0- The dashed red line marks 
Psus ~ 0.31%, the sustainable noise rate of the code, the lim¬ 
iting value of pth from to the fitting as TV — > oo. 

logical measurement. To perform the Z logical measure¬ 
ment [HISD, we measure each individual qubit of the 
code in the Pauli-Z basis. This destructive transversal 
measurement gives us the eigenvalues of stabilizers S q 
to diagnose E, and to thus recover the eigenvalue of Z. 

During readout, measurement errors and physical er¬ 
rors have an equivalent effect; both appear as bit flips. To 
simulate errors that occur during the readout process, we 
apply the noise operator E(p) to the encoded state before 
decoding. We therefore calculate logical failure rates 

Pfaii(A0 = prob (C(N)E(p)E £ G), (2) 

where C(N) = D(M 0 (E(p)\ipisr))). We evaluate Pf a ii(TV) 
values using Monte Carlo simulations. 

To analyze the performance of the proposed decoding 
scheme, we first look to find the sustainable error rate 
of the code, p sus , below which we can maintain quantum 
information for an arbitrary number of correction cycles. 
The discovery of such a point suggests that the error- 
correlations caused by our correction protocol do not ex¬ 
tend beyond a constant, finite, and decodable length, 
thus showing that we can preserve information indefi¬ 
nitely with arbitrarily high fidelity in the p < p sus regime. 

We define p sus as the threshold error rate, pth, at the 
TV —> oo limit, where the threshold is the error rate be¬ 
low which we can decrease -Pfail arbitrarily by increasing 
l ns]. In the inset of Fig. [4] we show the near thresh¬ 
old data we use to evaluate a threshold, where we show 
the data for TV = 8 as an example. We give the details 
of the fitting model we use to evaluate thresholds in the 
Methods Section of this Manuscript. 

We next study the evaluated threshold values as a func¬ 
tion of N. We show this data in the plot given in Fig. [5] 
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The data shows that p t h converges to p sus ~ 0.31% where 
we fit for values N < 8 . We obtain this value with a fit¬ 
ting that converges to p sus , namely 

Pth{N) = p sus [1 - (1 - pth(0)/psus) exp(- 7 lV)]. (3) 

We find pth(0) ~ 0.46% and 7 ~ 1.47. The convergent 
trend provides evidence that we achieve steady-state con¬ 
finement in the high-iV limit, as is required of a practical 
error-correction scheme. 

To verify further the threshold error rates we have de¬ 
termined, we next check that the logical failure rate de¬ 
cays as a function of system size in the regime where 
p < pth- hr the large N limit, we fit our data to the 
following hypothesis 

Pfaii(N) = {N + l)Aexp log (^~^j > ( 4 ) 

where A, a and /3 are positive constants to be determined 
and d is the distance of the code. We evaluate the vari¬ 
ables in Eqn. Q as A ~ 0.033, a ~ 0.516 and /3 ~ 0.822 
using ~ 5000 CPU hours with data for N < 10. Details 
on the fitting calculation are given in the Methods Sec¬ 
tion. We plot the fitted scaling hypothesis on Fig. [4] for 
the case of N = 8 to show the agreement of Eqn. Q with 
the available data. 

DISCUSSION 

To summarize, using only a simple decoding scheme, 
we have obtained threshold values that lie within an or¬ 
der of magnitude of the optimal threshold for the toric 
code under the phenomenological noise model. Moreover, 
we can expect that higher thresholds are achievable using 
more sophisticated decoding strategies (361 [37101 HU [521- 
[58] . It may be possible to achieve a sufficiently high sus¬ 
tainable noise rate to become of practical interest, thus 
meriting comparison with the intensively-studied surface 
code 59, 2>D]. To this end, further investigation is re¬ 
quired to learn its experimental viability. 

To continue such a comparative analysis, one should 
study the code using realistic noise models [23l that re¬ 
spect the underlying code hardware. We expect that the 
threshold will suffer relative to the surface code when 
compared using a circuit-based noise model where high- 
weight gauge measurements are more error prone [61] . 
Fortunately, gauge color code lattices are known where 
face operators have weight no greater than six [7]. While 
this is not as favorable as the weight-four stabilizer mea¬ 
surements of the surface code, given the ability to per¬ 
form single-shot error correction, and 7r/8-gates through 
gauge fixing, we argue that the gauge color code is de¬ 
served of further comparison. 

To the best of our knowledge, we have obtained the 
first threshold using single-shot error correction. Funda¬ 


mentally, our favorable threshold is achieved using redun¬ 
dant syndrome data to identify measurement errors. It is 
interesting to ask if we can make use of a more intelligent 
collection of measurement data to improve thresholds fur¬ 
ther. Discovering single-shot error-correction protocols 
with simpler codes might help to address such questions. 
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METHODS 
Threshold calculations 

The threshold error rate, p t h, is the physical error rate 
below which the logical failure rate of the code can be 
arbitrarily suppressed by increasing the code distance. 
We identify thresholds by plotting the logical failure rate 
as a function of physical error rate p for several different 
system sizes, and identify the value p = pth such that 
Pfaii is invariant under changes in system size. In the 
main text we show the data used for a specific threshold 
calculation in Fig. [4] where we use the N = 8 data as an 
example, and in Fig. [5] we evaluate threshold error rates 
for the gauge color code as a function of the number of 
error correction cycles, N. 

We evaluate threshold error rates by performing p ~ 
10 4 Monte Carlo simulations for each value of p close to 
the crossing point using codes of system sizes L = 23, 29 
and 35, except for the case that N = 0 where we evaluate 
the logical failure rate with system sizes L = 31, 39 and 
47. Simulating larger system sizes where N = 0 is possi¬ 
ble as in this case we read out information immediately 
after encoding it, such that E = 1 as shown in Eqn. (|2|. 
We therefore need not perform syndrome estimation in 
the N = 0 simulation. The threshold error rate at N = 0 
is thus the threshold error rate of the clustering decoder 
for the gauge color code where measurements are per¬ 
formed perfectly, i.e., q = 0. 
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We identify the crossing point by fitting our data to 
the following formula 

Pfail = B 0 + BiX + B 2 x 2 , (5) 

where x = (p — pthjL 1 ^ and Bj , p t h and p are constants 
to be determined. We show an example of this fitting in 
the Inset of Fig. [4j 

At the threshold error rate the code produces logical 
failures at a rate between 0.075 and 0.27 depending on TV. 
We expect such behavior as the number of logical failures 
will increase with repeated use of a decoder. We therefore 
obtain between ~ 750 and ~ 2700 logical failures per data 
point close to the threshold error rate. 

Overhead analysis 

Here we summarize the resource-scaling analysis we 
give in the p < p t h regime. We suppose the logical failure 
rate in this regime scales like 

-Pfaii(A0 « (TV + 1) A exp (a log(p/p t h (A0) d' 3 ), (6) 

where A, a and 0 are constants to be determined, d is 
the distance of the code, p is the error rate and TV is the 
number of uses of the decoder we make before readout. 
The value of p t h(TV) is determined by the method given 
in the previous Subsection. This equation is derived by 
assuming that a single use of the decoder will fail with 
probability Pf a ii = Aexp(ad^ log(p/p t h)) in the low-p 
regime. We then calculate to first order the probability 
that the decoder will fail a single time in TV + 1 uses to 
give Eqn. ([6]), where we include an additional use of the 
decoder to account for a possible logical failure during 
readout. 

We manipulate Eqn. © to show a method to evaluate 
a and 0. We first take the logarithm of both sides of 
Eqn. © to find the linear expression 

y = \og((N + 1)A) + ad p u, (7) 

where we write y = log Pf a n and u = log(p/pth)- We then 
take the gradient, g(d ) = dy/du, from Eqn. Q to find 

logff(d) = loga + /31ogd. (8) 

In Fig. [6] we plot logPf a ii as a function of logp. For 
fixed system sizes we observe a linear fitting that Eqn. Q 
predicts. We take the gradient of each of these fittings 
to estimate g(d), as given in Eqn. Q. Then, to find 
a and 0, we plot log g(d) as a function of logd where 
the gradients are taken from the fittings shown in Fig. [6] 
This data is shown in Fig. [7] We can now determine a 
and 0 using, respectively, the d = 1 intersection, and the 
gradient of a linear fit shown in Fig. [7] 

The logical failure rates we use to find a and (3 are 
found using between 10 5 and 10 6 Monte Carlo samples 



FIG. 6. Plot showing log Pf a n as a function of logp for the case 
that N = 8 where we have plotted system sizes L = 17, 23, 29 
and 35 shown in blue, yellow, green and red, respectively. The 
error bars show the standard error of the mean given by the 
expression A logp = [(1 - Pf a ii)Pf ai i/p] i/2 /Pf ai i where y is the 
number of Monte Carlo samples we collected. The logarithm 
of the gradients found for the linear fittings are plotted as a 
function of log d in Fig. [7] 

log(g) 



FIG. 7. Plot shows the available data fitted to the trend 
anticipated in Eqn. |8|. The Figure shows the logarithm of the 
gradients found in Fig. [6] log p(d), plotted as a function of the 
logarithm of the code distance, logd, for the case where N = 
8. Error bars show the standard error of the mean, evaluated 
using the LinearModelFit function in Mathematica. For 
this case we obtain a fitting log g(d) ss 0.81 logd — 0.63, as is 
shown in the plot. 

for each value of p where we only take values of p < 
0.8 x pth for each TV. We discard data points where we 
observe fewer than ten failures for a given p. The data 
was collected over 5000 CPU hours. 

We plot the values we find for a and /? as a function of 
TV in Fig. [8] Importantly, we observe convergence in the 
large TV limit. We see this using the fitting functions 

a(N) = a x [1 - (1 - a 0 /aoo) exp(- 7 a TV)], (9) 

and 


I3(N) = Poo [1 - (1 - 00/Poo) exp(- 7 /3 TV)], (10) 






a.iB 



FIG. 8. Unitless values, a and /3, plotted as a function of 
error-correction cycles, N, shown in blue and yellow respec¬ 
tively. The error bars show the standard error of the mean, 
which are found using the NonLinearModelFit function in¬ 
cluded in Mathematica. The a and /3 data points are fitted 
to Eqns. and ( |10| ) , respectively. The fitted functions are 
shown in the plot. 


to find the N —> oo behaviour of our protocol, where ao, 
A), c*oo, /3oo, 7 a and 7^3 are constants to be determined, 
such that 

Ooo = lim a(N), (3^ = lim (3(N). (11) 

AT—Kx> N—> 00 

We fit these functions to our data to find 

aoo = 0.516 ±0.005, A*, = 0.822 ± 0.004, (12) 

together with the following values ao = 0.65 ±0.02, 7 a = 
1.9T1.5, A) = 0.73±0.01 and 7 ^ = 1.7±1.4. The fittings 
are shown in Fig. [ 8 j 


A 



FIG. 9. Figure shows unitless values of A numerically deter¬ 
mined as a function of the number of error-correction cycles 
we perform, N. The error bars show the standard error of 
the mean which is calculated using Mathematica. We fit 
the data to Eqn. (|13|) as we show in the Figure. 


Finally, given that we have evaluated a(N) and j3(N) 
for different values of N, we use these values together 
with Eqn. (| 6 | to determine A as a function of N. We 
show the data in Fig. [9j We fit the values of A(N) to the 
following expression 

A(N) = Aoo [1 - (1 - Ao/Aoo) exp(- 7 A N )], (13) 

to find values A^ = 0.033 ± 0 . 001 , A 0 = 0.09 ± 0.01 
and 7 a = 2.1 ± 0.5, thus giving all of the variables, A ^, 
ttoo and Ax, we require to estimate the N —> 00 behav¬ 
ior of the decoding scheme in the below threshold regime. 
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SUPPLEMENTARY FIGURES 
Supplementary Figure 1 



(a) An odd and an even unit cube that each consist of five 
tetrahedra. The odd unit cube differs from the even unit 
cube by a ±7r/2 rotation about any of the canonical axes. 

(b) An exploded even unit cell that reveals the internal 
structure of a unit cube, (c) A fundamental tetrahedron 
of the lattice whose vertices are correctly four colored. 
We color each face of every tetrahedron with the color 
that is not the color of any of its three vertices. The face 
that is not visible is colored green. 


Supplementary Figure 2 



A cube of linear size L = 5 formed by stacking odd and 
even unit cubes. The odd and even unit cubes, shown in 
grey and white, are stacked such that no two even unit 
cubes meet at a face and no two odd unit cubes meet at 
a face. 


Supplementary Figure 3 



We remove tetrahedra from the cubic lattice shown in 
Supplementary Figure 2 to obtain the four-sided tetrahe¬ 
dral structure we require of the gauge color code lattice. 
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The figure shows some of the tetrahedra that have been 
removed from the cubic lattice to the right of the Figure. 
The remaining tetrahedra do not define a gauge color 
code dual lattice without replacing some of the removed 
tetrahedra, as we explain in Supplementary Note 1 and 
show in Supplementary Figure 4. 


The L = 5 lattice where the layers of tetrahedra are 
separated. Interior tetrahedra are colored white and grey. 


Supplementary Figure 6 


Supplementary Figure 4 




We modify the lattice of Supplementary Figure 3 to uni¬ 
formly color the boundaries by replacing some of the 
tetrahedra that were removed from the cubic lattice. 
Patches of six incorrectly colored faces are outlined with 
white hexagons. We add additional tetrahedra to these 
incorrectly colored patches such that all four boundaries 
are uniformly colored. The Figure shows some hexagons 
where the additional tetrahedra are already added. 


Supplementary Figure 5 



A repeating cell of the lattice that consists of eight of 
the fundamental unit cubes. Stabilizers of weight eight 
and thirty-two are represented by diamonds and circles 
on the vertices of the repeating cell respectively, as shown 
in the key at the right of the Figure. Opposite faces of 
the repeating cell are equivalent. 


Supplementary Figure 7 



Qubits on the exterior of the gauge color code lattice. 
Qubits lie on all the faces, on the edges where two differ¬ 
ently colored boundaries meet, fattened and colored in 
grey in the figure, and on the vertex where three differ¬ 
ently colored boundaries meet, marked by a large yellow 
vertex in the Figure. Qubits are numbered in accordance 
with the text in Supplementary Note 2. Qubits 7, 8, 13 
and 15 he on lattice edges and qubit 14 lies on the vertex 
where the red blue and green vertices meet. All other 
labeled qubits lie on the external faces of the lattice. 
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Supplementary Figure 8 



Errors introduced during fault-tolerant error correction, 
(a) An error string that, given perfect measurement out¬ 
comes, generates two blue stabilizer defects, (b) Mea¬ 
surement errors appear as strings on the dual lattice pic¬ 
ture, which are shown as multi-colored strings in the dia¬ 
gram. Gauge defects are identified where strings of mea¬ 
surement errors terminate. The original physical error 
is marked in grey, (c) During syndrome estimation, we 
incorrectly estimate the position of the measurement er¬ 
rors, and thus incorrectly identify the true position of 
the stabilizer defect. The correction operator we apply is 
represented by a blue string, (d) The correction operator 
we apply to correct the estimated stabilizer defect intro¬ 
duces a new error to the system. The net error that is in¬ 
troduced by the initial physical error and the inaccurate 
correction operator effectively acts like the discrepancy 
between the estimated position of the stabilizer defect 
from its true position. This discrepancy is marked by 
the blue error string in the diagram. 

SUPPLEMENTARY NOTES 
Supplementary Note 1 

Here we elaborate on the dual lattice we use to con¬ 
struct the gauge color code simulated in the main text. 
The gauge color code dual lattice is a neatly stacked 
structure of many tetrahedra. Conveniently though, the 
lattice dual to the lattice geometry shown in Figure 1(a) 
in the main text is composed of odd and even cubic units 
of five tetrahedra, shown in white and grey, respectively, 
in Supplementary Figure 1(a). We also show the de¬ 
composition of an even unit cell in Supplementary Fig¬ 
ure 1(b). An odd unit cell differs from an even unit cell 
only by a 7 t/ 2 rotation about any of the canonical axes. 
Unit cubes on the boundaries of the lattice are modified 
simply by removing subsets of their tetrahedra. 

To ease the explanation of the lattice construction, it 
is instructive to begin with a cubic block of fundamental 
cubes, as shown in Supplementary Figure 2. The block 
must have an odd linear dimension of unit cubes. Tetra¬ 


hedra are removed from this block to find an appropriate 
four-sided structure, and finally some of the tetrahedra 
that have been removed are once again replaced to find 
suitable boundaries that can be correctly colored. 

In the dual representation we require that the vertices 
of the lattice are four colorable, i.e. we can consistently 
color every vertex with one of four colors, red, green, yel¬ 
low or blue, such that no two vertices of the lattice that 
share an edge can have the same color. We show a single 
tetrahedra whose vertices are correctly colored in Sup¬ 
plementary Figure 1(c). The vertices of the dual lattice 
in Supplementary Figure 2 are indeed four colorable. 

The gauge color code requires four distinct boundaries 
where boundaries differ by the subset of colors of the 
vertices that lie on their surface. Specifically, a bound¬ 
ary of a given color contains no vertices of that color, 
i.e. a green boundary contains only red, blue and yellow 
vertices. We therefore look for a tetrahedral structure 
with four boundaries of four different colors. We design 
the correct structure by first removing large subsets of 
fundamental tetrahedra to find the global shape of the 
lattice. We find a four-sided tetrahedral structure by re¬ 
moving four corners of the cubic block of unit cubes, as 
shown in Supplementary Figure 3. 

The lattice shown in Supplementary Figure 3 does not 
define the gauge color code. This is because we cannot 
four-color the lattice such that all four boundaries of the 
lattice are correctly three colored, as we have explained. 
However, by replacing some of the tetrahedra we have 
removed from the lattice, we recover suitable boundaries. 
To find three-colored boundaries it is convenient to color 
the faces of the fundamental tetrahedra. We assign to 
each face of the tetrahedra the color that is not given to 
any of its adjacent vertices, as shown in Supplementary 
Figure 1(c). With this face coloring, we require that the 
lattice has four distinct, uniformly colored boundaries, as 
shown in the main text in Figure 1(b). 

We modify the tetrahedral lattice shown in Supplemen¬ 
tary Figure 3 to find a suitable lattice. In Supplementary 
Figure 4 we color the external faces of the fundamental 
tetrahedra to show the additional tetrahedra we must 
replace on the surface of the sheered lattice in Supple¬ 
mentary Figure 3. In particular, we observe hexagonal 
‘patches’ of faces which are inconsistently colored com¬ 
pared with the rest of the boundary. We outline some 
of these hexagonal patches in white. At the centre of 
each of these hexagons lies a single vertex of a color that 
is not suitable for the given boundary. To rectify this 
we replace all of the removed tetrahedra to the lattice 
that were originally touching the central vertex of each 
hexagonal patch. Upon doing so we recover the dual lat¬ 
tice shown in the main text. In Supplementary Figure 4 
we show some of the ‘bulbs’ of sixteen tetrahedra we rein¬ 
troduce to the lattice. Some of the bulbs have already 
been reattached in the Figure. For the convenience of 
the reader we show the lattice separated into layers in 
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Supplementary Figure 5. 

We finally remark that the lattice we consider is partic¬ 
ularly convenient for simulation because the stabilizers lie 
on the vertices of a cubic lattice. In Supplementary Fig¬ 
ure 6 we show a repeating unit of the lattice geometry on 
a cubic lattice. In the Figure stabilizers are represented 
by circular and diamond-shaped vertices of appropriate 
color as explained by the key. 

Supplementary Note 2 

Here we analyse the lattice by counting the total num¬ 
ber of qubits as a function of L. We also count the num¬ 
ber of stabilizers, and the number of gauge operators. We 
also discuss the qubit support of the stabilizers, and the 
gauge operators. The quantities evaluated in this Note 
serve as good sanity checks for readers that are repro¬ 
ducing the gauge color code lattice. 

Qubits 

The gauge color code is a complicated system where 
qubits are placed on tetrahedra, and on subsets of faces, 
edges and vertices of the dual lattice we described in Sup¬ 
plementary Note 1. A qubit is placed on each of the ex¬ 
terior vertices where three differently colored boundaries 
meet. We therefore have vertex qubits 

Qv(L)= 4. 

A qubit is also placed on exterior edges of the lattice 
where two differently colored boundaries meet. We have 

Qe{L) = 6 L, 

qubits on exterior edges of the lattice. 

The lattice has qubits placed on the exterior faces of 
the dual lattice. We find that there are 

Q f (L) = 9L 2 — 5, 

face qubits. 

Finally, the lattice has a qubit on each of its tetrahedra. 
We find 

Qt{L) = (5L 3 + 24L 2 -2 L- 24)/3, 
tetrahedron qubits. 

A gauge color code following our construction therefore 
has Q{L) = Q v (L)+Q e (L)+Qf(L)+Q t (L ) qubits, which 
is explicitly 

Q{L) = (5 L 3 + 51 L 2 + 16L - 27)/3, 


Stabilizers in the Dual Lattice 

Each vertex of the lattice supports two stabilizers. We 
find the number of vertices 

v(L) = ( L 3 + 12L 2 + 5 L- 6)/3. 

Each stabilizer represented by a given vertex acts on ev¬ 
ery vertex, edge, face or tetrahedron qubit that contains 
the respective stabilizer vertex. To make this statement 
more rigorous, we can specify each object on the lattice 
that supports a qubit by a list of vertices, v, that are 
contained by the respective object. A vertex qubit is 
specified by a single vertex, v = {v}, an edge qubit is 
specified by a pair of vertices, e = {iq, i^}, each triangu¬ 
lar face that supports a qubit is specified by three vertices 
/ = {vi,V 2 , U 3 }, and a tetrahedron is specified by its four 
vertices t = {tq, V 2 , V 3 , v±}. 

Given that the different objects of the lattice that sup¬ 
port qubits, known as simplices, can be uniformly de¬ 
noted by a lists of vertices of varying length, we are free 
to group all vertex qubits, edge qubits, face qubits and 
tetrahedron qubits into the set of qubits, < 2 , using this 
simplicial description. Written as simplices, vertex, edge, 
face and tetrahedron qubits only differ by the length of 
their list of vertices. Using this notation we can conve¬ 
niently write down stabilizer operators 

s *=n s v = n ( i4 ) 

QBv Q3v 

where we take the product of all qubits Q € Q that con¬ 
tain vertex v. We point out that we have used notation 
here that is inconsistent with the main text, as here we 
index stabilizers by vertices, v, whereas in the main text 
we chose the index c to represent cells that support sta¬ 
bilizers. This reflects the change from the primal to dual 
lattice notation. 

Gauge Operators in the Dual Lattice 

We finally find the number of face operators we use in 
the simulation in the main text, and explicitly describe 
the supports of the gauge operators on the dual lattice. 
Gauge operators are represented on edges, and on ex¬ 
terior vertices of the dual lattice. Instead of counting 
the number of gauge operators exactly, we count distinct 
supports of gauge operators. The number of supports 
are exactly the number of faces on the primal lattice. 
For each support, s, we have two gauge operators, Gf 
and Gf. On the dual lattice, gauge supports are uniquely 
represented by either an edge, e, or a pair (v, c) that con¬ 
tains a vertex, v, and a color c £ C = {r, g,y,b}, as is 
defined in the main text. 

To count the gauge supports, we first find the number 
of edges contained on the lattice, e(L). We find this 


qubits. 
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number using the Euler characteristic 


External vertex gauge operators are written 


X' 3 D =v{L)-e(L)+f(L)-t{L), (15) 


where % 3 D = 1 for a ‘ball-shaped’ triangulation, such 
as that which describes the gauge color code, and where 
v(L), f(L) and t{L) are the number of vertices, faces and 
tetrahedra of the lattice, respectively. We already have 
the number of tetrahedra t(L) = Q t (L), and v(L) is al¬ 
ready counted to find the number of stabilizers. We easily 
find the number of faces of the lattice f(L) using the fact 
that each interior face lies on the surface of two tetrahe¬ 
dra, and each tetrahedra has four faces. Following this, 
up to the exterior faces, we can regard each tetrahedra as 
contributing half to the face count per face of each tetra¬ 
hedra. We therefore obtain a contribution of ~ 2 Qt{L) 
faces for the tetrahedra of the lattice. To account for ex¬ 
terior faces, we must add an additional half-unit per exte¬ 
rior face of the lattice to count the total number of faces, 
giving the number of faces f(L) = 2 Q t (L) + Qf(L)/2. 
Explicitly, we have 


f(L) = (20 L 3 + 123 L 2 -8 L- lll)/6. 


Using f(L) and Eqn. (15) we find the number of edges 


e(L) = (4L 3 + 33L 2 + 2 L- 27)/2. 


Each edge of the lattice represents one gauge support 
that acts on all edge, face, and tetrahedron qubits that 
contain the respective edge. We use once again the sim- 
plicial notation to denote gauge supports represented by 
edges, e, such that 

Gf= n °e = n ( ie ) 

Q3{vi,v 2 } Q3{vi,v 2 } 

where we take the product over all qubits Q £ Q that 
contain both vertices of edge e = {^ 1 ,^ 2 }. We point out 
that no vertex qubit on the lattice can contain the two 
vertices of an edge, and therefore vertex qubits are not 
found in gauge operator supports associated to edges. 

Gauge supports are also represented by vertices v on 
the exterior of the lattice. Each exterior vertex represents 
either one, two or three gauge supports. For exterior 
vertices that represent multiple gauge operator supports, 
we specify the supports uniquely with a vertex, v, and a 
color c. 

Specifically, a gauge support of a given exterior vertex 
v contains a subset of the vertex, edge or face qubits that 
contain v. The subset of lattice objects Q 9 v are taken 
to be all of those that do not contain any vertices of one 
particular color, c. 

To make this statement rigorous using the simplex no¬ 
tation, we define the function C : v -9 C that returns 
the color of the vertex as defined by the choice of the 
four-coloring of the lattice. 


g ? v , c) = n n z q> ( i? ) 

Q3v \c Q3v\c 

where we use a shorthand notation l Q 9 v\c' to denote 
qubits Q that contain vertex v, but do not contain any 
vertices of color c, i.e. where C(v') c for all vertices 
v' £ Q. We point out that all tetrahedra contain one 
vertex of each color, and as such, no tetrahedra appear 
in any gauge supports represented by exterior vertices. 

In some cases there are exterior vertices with multiple 
nonempty subsets Q 9 v\c for different choices of c. In 
what follows we explicitly consider exterior vertices that 
contain one, two, and three gauge operator supports. We 
will see that exterior vertices in the middle of a bound¬ 
ary will contain only one gauge operator support, exte¬ 
rior vertices found where two different colored boundaries 
meet contain two distinct gauge operator supports, and 
vertices that lie where three different boundaries meet 
contain three different gauge operator supports. 

Exterior vertices in the middle of a boundary, far away 
from any edge or vertex qubits, denote one gauge oper¬ 
ator support. We explicitly consider the example of the 
gauge support associated to the blue vertex shown on the 
red boundary in Supplementary Figure 7. The gauge op¬ 
erator support associated to this vertex act on the qubits 
labeled 1, 2, 3 and 4. All qubits that contain the exterior 
vertex of interest contain a yellow, a green, and a blue 
vertex, but no red vertex. 

Vertices that lie at a point where two differently col¬ 
ored boundaries meet represent two distinct gauge op¬ 
erator supports. We consider the green vertex in Sup¬ 
plementary Figure 7 that lies where the blue boundary 
and the red boundary meet. One of the supports acts on 
the qubits labeled 3, 4, 5, 6, 7 and 8. Common to all of 
these qubits is that none of their faces or edges contain a 
red vertex. The other support acts on qubits 7, 8, 9, 10, 
11 and 12. Different from the first support we discussed 
that is associated to the green vertex, none of the faces 
or edges involved in this support contain a blue vertex. 

We finally consider the support of the gauge operators 
associated to the vertices that lie where three bound¬ 
aries meet, such as at the large yellow vertex shown in 
Supplementary Figure 7 that contains the qubit indexed 

14. This vertex denotes three supports, two of which are 
shown in the diagram. One support acts on qubits 6, 7, 
14 and 15, the qubits that contain the yellow vertex of 
interest, but do not contain a red vertex. The other sup¬ 
port acts on the qubits numbered 7, 12, 13 and 14, the 
lattice objects that contain the large yellow vertex, but 
do not contain any blue vertices. The third support as¬ 
sociated to the large yellow vertex acts on the green face 
qubit that contains the large yellow vertex which cannot 
be seen in the diagram, together with qubits 13, 14 and 

15. 
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We can now count the number of gauge supports we 
must measure to realize fault-tolerant error correction 
using the gauge color code. We count the number of 
exterior vertices v E (L) with the expression 

v e (L) = (9 L' 2 - l)/2. 

Then, using that the 6 {L — 1) vertices that lie where 
two boundaries meet represent two gauge operator sup¬ 
ports, and four vertices that lie where three bound¬ 
aries meet represent three gauge operator supports, we 
find the number of gauge operator supports G(L) = 
e{L) + v e (L) + 6 (L — 1) + 8, giving the number of gauge 
operator supports 

G(L) = 2 L 3 + 21 L 2 + 7L- 12. 

Twice this quantity gives the number of face operators 
we use to perform fault-tolerant error correction with the 
gauge color code to identify both Pauli-X and Pauli-Z 
type errors. 

Supplementary Note 3 

In the main text we identify the threshold error rate 
as a function of the number of rounds of error correction 
that are performed, N, to show that the threshold error 
rate is robust in the limit N —► oo. It is important to look 
at the single-shot error-correction scheme after repeated 
applications because, when one considers noisy measure¬ 
ments, the error-correction protocol will leave some resid¬ 
ual noise on the code. In general, this residual noise is 
not easily characterized. In fact, the noise introduced 
to the system is a complex function of the physical er¬ 
ror rate, the measurement error rate, and the choice of 
error-correction protocol. Given that the nature of the 
residual noise is not well understood, it is not clear that 
the error-correction protocol will be able to successfully 
deal with the residual noise after many cycles of error 
correction. 

In Supplementary Figure 8 we give an example of a 
mechanism that leads to the development of a correlated 
error. We show a physical error in Supplementary Fig¬ 
ure 8(a). The error is represented by a string. Given 
perfect measurements, we expect to observe two stabi¬ 
lizer defects at the left and right ends of the error string. 

For this example, some measurement errors occur in 
the vicinity of the right stabilizer defect when we at¬ 
tempted to learn the positions of stabilizer defects. The 
measurement errors, where face operators return incor¬ 
rect measurement outcomes, are represented by multi¬ 
colored strings in Supplementary Figure 8(b). Gauge 
defects lie at the points where the strings of measure¬ 
ment errors terminate. Given these measurement errors 


occur, we cannot be sure of the true location of the right 
stabilizer defect. Instead, we can make use of the posi¬ 
tions of the three gauge defects, and knowledge of the 
noise model, to attempt to determine the location of the 
stabilizer defect. 

Our ability to predict the positions of stabilizer de¬ 
fects depends on our choice of syndrome estimation al¬ 
gorithm. In Supplementary Figure 8(c) we show where 
the syndrome estimation algorithm mistakenly predicts 
measurement errors. The measurement errors that have 
been predicted are represented by strings that terminate 
at the gauge defects, and branch at the estimated stabi¬ 
lizer defect. The incorrect estimation made by the algo¬ 
rithm leads us to believe that the stabilizer defect is not 
in its true location, but is displaced to some estimated 
location. 

Next, we apply a correction operator to attempt to 
repair the initial physical error drawn in Supplementary 
Figure 8(a). However, using the estimated position of 
the right stabilizer defect, we apply a correction oper¬ 
ator that connects the left stabilizer defect to the esti¬ 
mated right stabilizer defect, as shown by a blue string in 
Supplementary Figure 8(c), thus introducing additional 
errors to the code. The initial error, and the applied 
correction operator is shown in grey in Supplementary 
Figure 8(d). 

Due to the topological nature of the string errors in the 
gauge color code, the net effect of the initial error and the 
correction operator equates to the discrepancy between 
the true position of the right stabilizer defect, and the 
position of the estimated stabilizer defect. The effective 
error is shown by the blue error string in Supplementary 
Figure 8(d). This is easily checked as we can continuously 
deform the grey string onto the blue string. In other 
words, in the gauge color code the grey error string is 
equivalent to the blue error string up to multiplication 
by gauge operators. 

In general it is not clear how able a syndrome estima¬ 
tion algorithm is for correctly predicting the positions of 
stabilizer defects. A bad syndrome estimation algorithm 
might displace stabilizer defects over long distances com¬ 
pared to their true positions, in which case large corre¬ 
lated errors that we cannot correct may develop. More¬ 
over, it is not clear that the character of the residual 
noise remains constant over many error-correction cycles. 
Indeed, one should be concerned that an error correc¬ 
tion protocol might cause correlations to develop over 
repeated use of the error-correction procedure while in¬ 
formation is stored. It is therefore important to study 
a single-shot error-correction protocol over many cycles 
of error correction to interrogate its performance, and to 
check that the noise incident to a code achieves a steady 
state in the long-time limit. In doing so, we are able 
to establish that a code has a finite threshold after an 
arbitrarily long time. 



